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The effects of the electron-electron interactions in a graphene layer are investigated. It is shown that 
short range couplings are irrelevant, and scale towards zero at low energies, due to the vanishing of 
density of states at the Fermi level. Topological disorder enhances the density of states, and can lead 
to instabilities. In the presence of sufficiently strong repulsive interactions, p-wave superconductivity 
can emerge. 
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I. INTRODUCTION. 

Recent experiments pj|| report the existence of fer- 
romagnetic and superconducting fluctuations in graphite 
at unexpectedly high temperatures (T ~ 100 — 300K). 
The coexistence of both types of fluctuations suggests a 
common electronic origin for them. 

Motivated by these observations, we present here a 
study of the possible electronic instabilities of a single 
graphene sheet. Isolated graphene has the convenient 
property that the electronic states near the Fermi level 
can be described in simple terms. By symmetry, the 
lower and upper bands touch at the corners of the hexag- 
onal Brillouin zone. Near these points, the dispersion 
relation is isotropic and linear, eg = ufI^I, where vf is 
the Fermi velocity. The density of states at the Fermi 
level is strictly zero, and it rises linearly in energy. An 
effective long wavelength description of these electronic 
states can be written in terms of the Dirac equation in 
two dimensions (see below). 

The fact that a single graphene sheet is a semimetal 
modifies significantly the screening of the Coulomb in- 
teraction H. An effective low energy hamiltonian can 
be written, which can be treated by Renormalization 
Group methods It can be shown rigourously that 

the Coulomb interaction is a marginal interaction, which 
scales to zero at low energies or long wavelengths. At 
intermediate scales, however, the quasiparticle lifetime 
does not follow the usual e 2 dependence of Landau's the- 
ory of a Fermi liquid, but scales as |e| |6|, in agreement 
with experiments [Q. The RG approach is, in principle, 
valid in the weak coupling regime, (e 2 / (eo^f) "C 1, where 
e is the electric charge and eo is the dielectric constant. 
By using a RPA summation of diagrams, it can be shown 
that the low energy properties are not changed through- 
out the entire range of couplings [|| . 

The previous work mentioned earlier analyzed the 
the small momentum scattering due to the long range 
Coulomb interaction, as it is the only one which leads to 
logarithmically divergent perturbative corrections. Some 
electronic instabilities, like anisotropic superconductiv- 
ity, requires the existence of short range interactions with 
significant strength at finite wavevectors. We analyze 



in this work the role of these interactions in inducing 
instabilities of the electronic system. The next section 
describes the model. Then, the Renormalization Group 
equations for the different iteractions are written. In sec- 
tion IV, the role of topological disorder is analyzed, as it 
can lead to changes in the density of states which modify 
the scaling equations obtained earlier. The main conclu- 
sions are presented in section V. 

II. THE MODEL. 
A. Intralayer couplings. 

We analyze the low energy properties of a graphene 
sheet. We will only consider the modifications due to in- 
teractions and disorder in the low energy properties of the 
system. Thus, we need to describe the low energy elec- 
tronic states. A graphene sheet has an hexagonal sym- 
metry with two atoms per unit cell. The carbon atoms 
have four valence orbitals. Three of them build the sp 2 
bonds which give rigidity to the structure. The third 
orbital gives rise to the valence and conduction bands. 
These bands touch at the two inequivalent corners of the 
Brillouin zone (see Fig.JlJ]). From symmetry considera- 
tions, these bands are isotropic, and depend linearly on 
the wavevector. 

It can be shown that, in the long wavelength limit, the 
electronic wavefunctions near the corners of the Brillouin 
Zone are well described in terms of the two dimensional 
Dirac equation. Each of the two inequivalent points re- 
quires two Dirac spinors, each of them with its spin index. 
In the long wavelength limit, the The Fermi velocity, vf, 
can be expressed in terms of the matrix elements between 
nearest neighbor tt orbitals, t, as vf = (3ta)/2, where a 
is the C-C distance. 

Because of the collapse of the Fermi surface to isolated 
points, the kinematics are much simpler than the corre- 
sponding analysis for two "hot spots " in a square lattice 
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The hamiltonian is: 
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where a x and a y are 2x2 Pauli matrices. We have 
separated the long wavelength part of the Coulomb in- 
teraction from other possible short range interactions. 
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FIG. 1. Sketch of the Brillouin Zone of a graphene sheet, 
and the band dispersion near the Fermi energy. 



The couplings, gis'-,s,s' can be classified in an analo- 
gous way as in one dimension. The possible scattering 
processes are shown in Pig.g]. 

Because of the linear dispersion of the electronic states, 
we can use vp to transform time scales into length scales. 
Then, we can express the dimensions of all physical quan- 
tities in terms of lengths. Within this convention, we find 
that the dimension of the electronic fields is [iff] = l^ 1 , 
where / defines a length. A naive power counting anal- 
ysis shows that the Coulomb potential defines a dimen- 
sionless, marginal coupling, while the g's scale as I, and 
are irrelevant at low energies. This effect can be traced 
back to the vanishing density of states at the Fermi level. 
When a single Hubbard intrasite repulsion U is consid- 
ered, all interactions between electrons of opposite spin 
in eq. {[j]} are equal to lift, where f2 is the area of the 
unit cell, and the interactions between electrons of par- 
allel spin are zero. 
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FIG. 2. Possible short range couplings between electrons 
near the two Fermi points in a graphene layer. The full and 
broken lines correspond to electrons in the vicinity of each of 
the two Fermi points, a) Intra-singularity scattering (gintra). 
b) Inter-singularity scattering (<?i n tcr). c) Exchange scattering 

(^exchange ) ■ 



B. Interlayer couplings. 

So far, we have restricted our analysis to processes 
within an isolated graphene sheet. Neighboring layers 
are always coupled by the Coulomb interaction. In the 
following, we will neglect interlayer hopping, so as to 
be able to describe the electronic levels in terms of the 
Dirac equation, but we include the effects of the long 
range Coulomb interactions between layers. The inter- 
layer couplings give rise to the screening of the bare in- 
tralayer electron-electron interaction. We will treat these 
effects within the RPA, as depicted in Fig.||, following 
the analysis in M. 
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The intralayer interaction becomes: 
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where d is the distance between layers, and xo is the 
electron susceptibility of a single layer, given by: 
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FIG. 3. Diagrammatic representation of the Random 
Phase Approximation applied to the interlayer Coulomb in- 
teraction. The diagrams show the screening of the interaction 
between two electrons in layer I due to the polarization of lay- 
ers I', I".... 

The interlayer interactions are only effective when 
|q|d <C 1. Hence, if the lattice constant a is such that 
a <C d, they do not affect significantly the couplings be- 
tween electronic states in different Fermi points. 



III. SCALING ANALYSIS. 

In II it was shown that the electrostatic coupling, de- 
fined as e 2 /(eQVF), scales towards zero at low energies, for 
all values of the interaction. On the other hand, the ex- 
istence of scattering processes between the two inequiva- 
lent Fermi points can lead to instabilities at intermediate 
couplings. Different combinations of couplings lead to 
each instability. The system becomes ferromagnetic for 

Sufficiently large values of S'intra_L+fl , i I iter_L-fl , intra|| -flinter|| i 

where the subscripts || and _!_ denote the relative orien- 
tation between spins. An antiferromagnetic instability is 
driven by g in t ra ± + ffmterj. - <7mtra|| - ffintcr|| • The super- 
conducting phases can be s and p wave, depending on the 
relative phase of the gap at the two inequivalent points. 
However, for each k near the Fermi points, there are two 
electronic states, so that an additional index can be de- 
fined in the superconducting order parameter. Writing 
these two states as a two component spinor, we can write, 
in general: 



A,T,k 
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where a and b are constants. When the interaction is re- 
pulsive, the p-wave symmetry is favored (Ag = — A_g), 
as in a two dimensional electron system with two inequiv- 
alent van-Hove singularities at the Fermi level ||. The 
corresponding coupling is g intC r± + <?i n terJ_ - ^exchanged. - 
5exchango± ■ The diagrams which define the flow of these 
couplings are depicted in Fig.Q. 
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FIG. 4. Renormalization Group equations in the Cooper 
channel with p-wave symmetry, a) Sketch of the order pa- 
rameter in the Brillouin Zone, b) Diagrams involved in the 
calculation. 



The corresponding equations for the dimensionless ver- 
tices, T's, can be written as: 



cTi, 
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where we are omitting spin and flavor indices for simplic- 
ity, and dp is the (anomalous) dimension of the vortex, 
which includes, among others, the effects of the wave- 
function renormalization of the fields. To lowest order, 
df = 1. The first term in the r. h. s. of eqs.|5| is lin- 
ear, and it is absent in the flow of the couplings in the 
Cooper channel in a conventional metal. It reflects the 
irrelevance of these couplings in a semimetal. 

The flow in this channel becomes relevant if r cxc h a ngc > 
Tintcr and the values of the T's are of order unity. Note 
that the cutoff is assumed to be A w vp/a, where a 
is a length of the order of the lattice constant. The 
dimensionful inter-Fermi points and exchange couplings 



induced by the Coulomb interactions are <?j 



7M. 



Hence, the bare vortices, f q ~ e 2 /(eaVp). For reasonable 
values of eo ~ 4 — 8, this combination is, indeed, of order 
unity. 



IV. INFLUENCE OF DISORDER. 
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A. Topological disorder. 

The formation of pentagons and heptagons in the lat- 
tice, without affecting the threefold coordination of the 
carbon atoms, lead to the warping of the graphene sheets, 
and are responsible for the formation of curved fullerenes, 
like Cgo- They can be viewed as disclinations in the lat- 
tice, and, when circling one such defect, the two sub- 
lattices in the honeycomb structure are exchanged (see 
Fig.§). 




FIG. 5. Formation of a pentagonal ring in the honeycomb 
lattice. Points a,b,c,d... have to be identified with points 
a , b' , c', d! .... The defect can be seen as a disclination, defined 
by the straight dashed lines. 

The two fermion flavors defined in eq.{|I]} are also ex- 
changed when moving around such a defect. The scheme 
to incorporate this change in a continuum description was 
discussed in [jlO| . The process can be described by means 
of a non Abelian gauge field, which rotates the spinors in 
flavor space. The vector potential is that of a vortex at 
the position of the defect, and the flux is ±7r/2. 

Dislocations can be analyzed in terms of bound discli- 
nations, that is, a pentagon and an heptagon located at 
short distances, which define the Burgers vector of the 
dislocation. Thus, the effect of a dislocation on the elec- 
tronic levels of a graphene sheet is analogous to that of 
the vector potential arising from a vortex-antivortex pair. 
We can extend this description [jll), and assume that a 
lattice distortion which rotates the lattice axis can be 
parametrized by the angle of rotation, 9(r), of the local 
axes with respect to a fixed reference frame. Then, this 
distortion induces a gauge field such that: 



A(r) = 3V0(r) 



-i 

1 



(6) 



Thus, a random distribution of topological defects can 
be described by a (non abelian) random gauge field. The 
nature of the electronic states derived from the two di- 
mensional Dirac equation in the presence of a gauge field 
with gaussian randomness has received a great deal of 
attention, as it also describes the effects of disorder in 
integer quantum Hall transitions |l2[. The disorder is 



defined by a single dimensionless quantity, A, which is 
proportional to the average fluctuations of the field: 



(A(r)A(f)) = Ac5 2 (r-f) 



(7) 



It is known that A gives rise to a marginal perturbation, 
which modifies the dimensions of the fermion fields and 
enhances the density of states at low energies. A variety 
of analytical (l5j and numerical techniques jl6) has been 
used to study this problem. We will follow the Renor- 
malization Group scheme presented in fll2| . 

We first analyze the statistical properties of the gauge 
field induced by topological defects. Let us assume that 
the graphene sheet is warped, and that there is a random 
distribution of pentagons and heptagons, with density no 

— 1/2 

and average distance equal to Iq = n Q . The fluctua- 
tions in the gauge field induced by this distribution at 
a given point can be calculated by considering the effect 
of all defects located at distances between r and r + dr 
(see Fig.||), where r ^> Iq. The number of defects of 
each type is 2nrdrnQ. The angle, </>, which defines their 
position is a random variable. The contribution to the 
x-component of the gauge field from these vortices is: 



4(r = 0) 



r 

$ V 1 



(8) 




FIG. 6. Effect of vortices located at distances between r 
and r + dr from the origin (see text for discussion). 

where $o is the flux associated to a single vortex, and 
there is a similar equation for A y (r = 0). We now must 
integrate this value from Iq to R, where R is the radius 
of the sample. We obtain: 



4 



|A(r = 0)| 2 = 2vrno$olog 



(9) 



We can assume that the vector potential at positions sep- 
arated by distances greater than Iq are not correlated. 
Then, from eq.{Q}, we find: 



A = 2;r<I>n !<>;.{ ^ 



(10) 



which diverges slowly with the size of the system. The 
previous estimate assumed that the layers had a signif- 
icant amount of curvature at distances smaller than Iq. 
We can alternatively assume that pentagons and hep- 
tagons are bound in dislocations with average distance 
b. The vector field of a vortex-antivortex dipole decays 
as r~ 2 . A similar analysis to the one leading to eq.{|10|} 
gives: 



A oc ^Indisib 2 



(11) 



where ridisi is the density of dislocations. 

We will now assume that random fields induced by 
topological defects have the same statistical properties 
to those with gaussian disorder with the same value of 
A, which is the second moment of the distribution in 
both cases. Then, we can perform the Renormalization 
Group analysis discussed in |l2] ]. To lowest order, we find 
an interaction between fcrmion fields in different replicas 
of the type: 

S int = AY^ I [*A(?,*l)*B(?,*l)] m 

x [y B (r,t 2 )y A (r,t 2 )] n dt 1 dt 2 dr (12) 

where m and n are replica indices. This interaction leads 
to a logarithmically divergent self energy, which can be 
interpreted as a renormalization of the density of states 
p2| . We can include the corrections induced by the self 
energy in a renormalization of the wave function, giving 
rise to a change in the scaling dimension of the fields: 



2d* 



1 = 1-* 

7T 



(13) 



This expression has to be inserted in eq.{[s]}, modifying 
the flow of the couplings. 

The same result can be reached by analyzing the self 
energy corrections using standard techniques in the study 
of disordered electrons in arbitrary dimensions ]l3|] . To 
lowest order, the first correction to the Green's function 
is shown in Fig.0. This diagram leads to a self energy: 



E(F,P,w) « (G (r-r', w )A(r)A(P)) + 
= AG (r - ? , lu)S 2 (r - r') + • 



(14) 



where Go is the unperturbed Green's function. The real 
part of Go behaves as Go ~ to log(A/cj). Finally: 



2d* - 1 = 1 - — r I — I 
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FIG. 7. Diagrams contributing to the Green's function in 
a disordered electron system. 

The previous perturbative analysis can be generalized 
to arbitrary couplings, by mapping the non interactiong 
fermion problem in two spatial dimensions onto an inter- 
acting problem in 1+1 dimensions 0]. At energies below 
a scale A ~ Aexp[— 7r/(2A)], the backscattering between 
the two Fermi points leads to a scaling dimension which 
is independent of the disorder, d* = 1/7. 



B. Substitutional and site disorder. 

It can be shown that substitutional and site disorder 
can be incorporated into the Dirac equation through a 
change in the local chemical potential and the appear- 
ance of a mass term [lq| . Disorder of these types, with a 
gaussian distribution, defines a marginally relevant per- 
turbation [ji5| . Within the perturbative RG scheme de- 
scribed in the previous subsection, it can be shown that 
this perturbation leads to logarithmic corrections to the 
site-diagonal self energy, which can be incorporated into 
a renormalization of the wavefunction. Moreover, devia- 
tions from a gaussian distribution are relevant perturba- 
tions, which modify the results. 

Thus, while to lowest order substitutional disorder 
shows similar characteristics as topological disorder, 
higher order corrections lead to significant modifications, 
which, in addition, depend on the type of disorder Jl7| . 
For instance, a sharp divergence of the density of states 
has been found pq| , or a suppression of the density of 
states at low energies Ji~9[ | . 

In graphite, we do not expect a high concentration of 
charged impurities, which will lead to a strong site disor- 
der. Randomness in bond lengths leads to non diagonal 
disorder, which can be included in the topological dis- 
order discussed in the previous subsection. Hence, we 
expect diagonal disorder in pure graphite to be small, 
leading to minor corrections to the dependence of the 
density of states on energy, even when non perturbative 
terms are included Jl7| . 
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V. DISCUSSION. 
A. Other effects not included in the model. 

We have neglected the effects due to phonons, al- 
though, in principle, they can be incorporated into the 
framework used here. Our main purpose is the study 
of instabilities towards ground states which exhibit mag- 
netism or anisotropic superconductivity. We assume that 
the electron-phonon interaction will not change qualita- 
tively the possible existence of these instabilities. 

We have also not consider the interlayer hopping, t±. 
Within the RG scheme, coherent interlayer hopping is a 
relevant perturbation |20(| , leading to three dimensional 
behavior at low energies or temperatures. On the other 
hand, due to the vanishing of the density of states of 
a graphene layer, incoherent hopping between layers is 
irrelevant (note that it is a marginal perturbation in sys- 
tems with a finite density of states |2l| ) . 

In the presence of coherent interlayer hopping, our 
analysis is valid only at scales higher than t±, which 
has been estimated, by band structure calculations to be 
t± w 0.27eV p^|. This bare value will be reduced by the 
many body effects, and the wave function renormaliza- 
tion considered here. However, in a perfect system, the 
validity of our calculations are limited to a range between 
t « 2.4eV and the renormalized value of t±. 

The coherent interlayer hopping can modify our results 
in various ways: i) The coupling between layers induces 
a crossover to 3D behavior, enhancing the 2D instabili- 
ties discussed here, ii) The dispersion of the electronic 
bands in the third dimension leads to the existence of 
small electron and hole pockets, increasing the density 
of states. If the couplings are not modified, this finite 
density of states will also strenghten the instabilities, ii) 
The density of states at the Fermi level induces metallic 
screening, and changes the interactions at low energies. 
It is unclear to us how our results are modified in case 
iii)- 

Our calculations have a wider range of validity in the 
presence of disorder, where coherent hopping over dis- 
tances longer than the electronic mean free path is sup- 
pressed. As mentioned earlier, incoherent local hop- 
ping can be considered an irrelevant perturbation which 
should not modify qualitatively the results presented 
here. 



B. Analysis of the couplings. 

Our analysis considers the role of electron-electron in- 
teractions in a graphene layer. Spin dependent inter- 
actions, like a Hubbard on site term, naturally lead to 
magnetic phases. In the absence of disorder, a minimum 
value for the Hubbard repulsion is required before the 
onset of antiferromagnetism p3|, in agreement with the 



analysis presented here. This phase, however, lacks ex- 
perimental confirmation. It is also known that, within 
the Hartree-Fock approximation, a nearest neighbor re- 
pulsion, V, induces a charge density wave ground state, if 
U-W < 0, and U and V are sufficiently large In 
addition, there is a region in the phase diagram where U 
and V almost cancel, leading to a paramagnetic ground 
state. Realistic values of these parameters sugest that a 
graphene layer lies in this region p4|-p(| . It is reasonable 
that longer range correlations can make this state un- 
stable. These calculations do not consider longer range 
interactions. For decoupled graphene layers, the vanish- 
ing of the density of states at the Fermi level leads to the 
absence of metallic screening, so that spin independent, 
long range interactions are expected. 



C. Low temperature phases. 

We have considered the possibility of ferro- and an- 
tiferomagnetism, and p wave superconductivity as the 
most likely low temperature phases. The competition 
between them depends on the spin dependence of the in- 
teractions. Spin independent couplings favor supercon- 
ductivity while a strong spin dependence, like the on 
site repulsion of the Hubbard model, will lead to a mag- 
netic ground state. Finally, the relative stability of ferro- 
and antiferromagnetism depends, among other things, on 
the existence of an underlying bipartite lattice. In the 
presence of a sufficiently strong topological disorder, we 
expect that ferromagnetism will prevail over antiferro- 
magnetism, as the existence of pentagons and heptagons 
leads to the frustration of antiferromagnetic order. The 
same argument can be applied to the charge density wave 
state considered in p4p5| . 

A detailed study of the competition between ferro- 
magnetism and p-wave superconductivity lies beyond the 
scope of this work. It depends on the balance between 
the on site, spin dependent interactions, and the longer 
range, spin independent couplings. Ferromagnetism is 
favored by the existence of a sufficiently strong forward 
scattering between electrons of opposite spin, at momen- 
tum transfer q w 0. This coupling depends on the nature 
of the screening, which, in turn, depends on the density 
of states near the Fermi level, and on the degree of disor- 
der. On the other hand, if the main interactions are spin 
independent, ferromagnetism will be suppressed, and the 
leading instability is p-wave superconductivity. 



D. General features of the possible superconducting 
instability. 

In the following, we will focus discuss some qualitative 
features of the superconducting transition. A quatitative 
estimate of the critical temperature is beyond the scope 
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of our RG scheme, although we can discuss the depen- 
dence of T c on various quantities. 

It is first interesting to note that superconductivity at 
low temperatures was observed in graphite intercalation 
compounds p7[ . The origin of this superconductivity 
is not completely understood. The critical field shows 
an anomalous dependence on temperature pq , unlike in 
conventional s-wave superconductors. This dependence 
has been explained in terms of a two band model f p9| . 
This model is similar to the two point model discussed 
here, except that the two bands considered in |2^] corre- 
spond to a carbon and a dopand band. The temperature 
dependence of the critical field should be, however, sim- 
ilar in the two cases. 

The critical temperature at which an instability de- 
scribed by eq.{|]} sets in is: 



T r = A 



r 



r 



(16) 



where T is the appropiate vortex required to drive the 
instability. There is a transition if T > T c = df.. For 
df = 0, this expression reduces to the usual BCS for- 
mula, T c = Aexp(— 1/fo), and f c = 0. The disorder 
influences the scaling of the fermion fields, d<s,, which, in 
turn, modify df,: 



df = Ad^ — 3 = 1 — 



2A 



(17) 



where A is given in eq.{10} or eq.{|Tl]}. The critical tem- 
perature depends exponentially on the disorder. The ex- 
pression in eq.{|l6|} is only valid if T c -C vp/d, where d is 
the typical distance above which eq.{10} or eq.{|il"|} hold. 

We can make a simple estimate of the role of disorder 
by assuming that, for certain average separation between 
defects, Iq, dp = 0, and the value of the critical tempera- 
ture is T c max . Then, if the disorder is reduced, we expand 
on dg -C 1, and we obtain: 



the optimal concentration of defects will be in the range 
lo ~ £o- Assuming that T™ ax ~ 300K, this estimate gives 
for the mean distance between defects Iq ~ 30 — lOOA. 



E. Origin of disorder in graphene sheets. 

It is known that electronic properties of graphite, like 
the resistivity, are sample dependent p0[ , and localiza- 
tion effects due to disorder have been observed As 
discussed in section IV, the effect of topological disor- 
der depends on whether the graphene sheets present a 
finite density of disclinations, leading to corrugated and 
warped surfaces, or the main source of disorder is due 
to dislocations. Aggragations of graphite nanoparticles 
of polyhedron shapes, whose curvature is not completely 
characterized are discussed in Warped layers, with 
curved regions which are reminiscent of the spherical 
fullerenes have been observed |33| . These structures seem 
similar to proposed models of negatively curved graphene 
layers |3^] . Theoretically, these compounds (schwarzites) 
are supposed to be very stable, and contain a macroscopic 
fraction of heptagonal rings. A material with these char- 
acteristics is probably best described by a random distri- 
bution of disclinations, with mean separation equal to a 
few lattice spacings. Calculations of the electronic den- 
sity of states of the model proposed in |34| show that 
it loses the semimetallic properties of graphite, in agree- 
ment with the discussion here ]35| . Compounds with 
these characteristics can be good candidates for intrinsic 
p-wave superconductivity. 

From the difference between eqs.{|l( 



and {11}, it is 
clear that a non corrugated graphene sheet has a much 
lower density of states than a significantly warped one, 
and a reduced tendency towards electronic instabilities. 
In highly disordered graphite, however, it is possible that 
regions with different degrees of corrugation coexist giv- 
ing rise to the behavior reported in . 



T c u « Ae 

^ rjvm 



(18) 



where k 



oc 



is a numerical constant and I is the aver- 



age distance between defects. The superscript stands 
for the fact that frustration effects in the superconducting 
phase are not considered (see below) . Finally, we can get 
a rough estimate for Iq by considering that a sufficiently 
large concentration of defects leads to pair breaking and 
reduces T c in an anisotropic superconductor. The reduc- 
tion of T c is given, approximately, by |1 1| : 



(19) 



where £o = v f/T® is the coherence length of the super- 
conductor, and c is a constant of order unity. Hence, 
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